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Abstract 

Quantum mechanical scalar particle with polarizability is considered in the presence of the 
Coulomb field. Separation of variables is performed with the use of Wigner /^-functions, the 
radial system of 15 equations is reduced to a single second order differential equation, which 
among the Coulomb term includes an additional interaction term of the form era 2 /M 2 r . 
Various physical regimes exist that is demonstrated by examining the behavior of the curves 
of generalized squared radial momentum P 2 (r). Eigenstates of the equations can be con- 
structed in terms of double confluent Heun functions. Numerical analysis proves the existence 
of the bound states in the system; the lowest energy level and corresponding solution are 
calculated based on generalization of Ritz variational procedure. 



1 Separation of the variables 

A particle with spin and polarizability can be described with the use of tetrad formalism by 
the following equation (the main references to original papers concerning this model are given 
in ID) 

[T a (x) (d a + Bjx)) -m]*(x) = 0, 

1 



T a (x) = T a e? a) (x) , B a (x) = ^ J ab e?V a (e m ) . 



(a) 



(1.1a) 



In Minkowski space and in spherical coordinates and tetrad, eq. (1.1a) takes the form 

9(x) = Q, (1.16) 
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(1.1c) 



General form for 15-component spherical wave function e,j,m is (more details see in [2]; the 
notation for Wigner functions is used; D a = D J _ m a ((fr,6,0)): 



C(x) = e~ tet C(r)D , C (x) = e~ iet C (r)D , * (a?) = e^o^A) , 



C{x) 



-icf 



Ci(r) D_! 
C 2 (3) D 
G 3 (r) D +1 



<&(x) 



-let 



*i(r) D-i 
$ 2 (r) D 
Mr) D +1 
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E(x) 



-let 



Ei(r) D_! 
E 2 {r) D 
E 3 (r) D +1 



, H(x) 



-let 



Hi(r) £>_i 
H 2 (r) D 
H 3 (r) D +1 



After separation of variables we obtain the radial system [2] 



-i(e+-)C - + -)C 2 - -(d + C 3 ) = mC, 
r dr r r 

-i{e+-)C - «r(-£- + ~)E 2 - — (E 1 + E 3 ) = mC 
r dr r r 

i{e + -)(j Ei + o- ( — + -)H 1 C + H 2 = mCi 

r ar r r r 



• / a\ ^ d _ 
i (e+ -) a £ 2 + -r C 
r ar 



r 

n , , d 1 s „ £/ „ i f cr 



i (e + -) a E 3 - a( — + - ) H 3 - - C 
r ar r r 



(Fx - H 3 )=mC 2 , 

H 2 = mC 3 



ot > 



-i (e + — ) C = m $ 
r 



- (7 = m <3?i 

r 
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C = m $2 , - 


— C = m 
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(±)[- 


-» (e+ -) $i + 
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m £?i , 


(±)[- 


i (e+ -) $ 2 - 
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aV $o] = 


= m E 2 , 


(±)[- 


-i (e+~) $ 3 + 
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^*o] = 
r 


m £?3 , 


±)[-i 


/ d 1 . , 
ar r 


$ 2 

r 


= m i?i 



(±) - ( $a - $ 3 ) = m F 2 , 
r 

(±) [ +i ( + - ) $ 3 + — <*>2 ] = m H 3 . 

ar r r 

From (1.2c) it follows that $ 3 = Then from (1.2d) we get E 3 = +E U H 3 

Finally, from (1.26) it follows C 3 = + C\. So, the collected restrictions are 

$ 3 = + *i, C 3 = + Ci, E 3 = + Ex, H 3 = -Ht, H 2 = 

Taking into account (1.3), the radial system reads 

-i (e + -) Co - ( 4- + -) C 2 - - 2Ct = m C , 
r dr r r 

-i (e + -) C - a ( 4- + - ) ^2 - — 2E 1 = mC , 
r dr r r 

i(e + -)a Ei + a ( — + - ) iTi - - C = m Ci , 

r ar r r 



i (e + -) a E 2 + 4- C - — 2H 1 =mC 2 , (1.46) 
r dr r 

-i (e + — ) C = m $o , C = m $i , 

r r 

A C = m $ 2 , (1.4c) 

dr 

cx d 
(±) [ -i (e + -) $ 2 - -J- $o ] = m £ 2 , 

r or 

(±) [ -t ( + - ) $i - — $ 2 =mH l . (lAd) 
dr r r 

Taking expressions for Cj according to (1.4c) , from (1.4d) one gets 

E 1 = 0, E 2 = (±)(-A^)C, H 1 = 0. (1.5) 

Now, from (1.46) with the use of (1.5), one obtains 

_ .. a. , iacr dC v dC , acr .., „. 

mC = -i(e + -) ± — 5-- "r" > mC i = C . mC 2 = -r ± — 5-~C . (1.6) 

r m^r ar r ar m z r z 

Finally, we arrive at a second order differential equation C(r) (changing to 2 to— M 2 , and 2z/ 2 to 

i(i + i)): 

ydr 2 r dr r r 2 M 2 r A J 

2 Qualitative analysis of the radial equation 

With the use of notation (the signs ± can be included into the parameter a) 

e 2_M 2 = -K 2 , j(j + 1) - a 2 = J 2 , a ^p =Y? (2- 1 ) 
the main equation (1.7) reads 

C(r)-lm, g + (-^2= _£ + *),„)_„. ,2,, 

Let us examine the behavior of the squared radial momentum. Near the origin and at infinity 
we have 

P 2 (r^0) r , P 2 (r -»• oo) ~ (e 2 - M 2 ) . 

r 4 

To describe classical turning points it is convenient to factorize the expression for P 2 (r) 

2 , , (e 2 - M 2 )r 4 + 2ear 3 - J 2 r 2 + S 2 

P V) = ~4 = 

= (e 2 - M 2 )(r - ri)(r - r 2 )(r - r 3 )(r - r 4 ) = Q ^ ^ 
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Figure 1: Finite classical motion: r £ [r 3 ,r 4 ] 



The roots of the 4-th order polynomial obey relations below 

2ea 



J 2 



e 2 _ M 2 ^ r l + r 2 + r 3 + n , 

riV2 + nr 3 + r±V4 + r 2 r 3 + r 2 r 4 + r 3 r 4 , 



(2.4a) 

(2-46) 
(2.4c) 

e 2 _ y[2 — • ■ (2 Ad) 

First, let us consider the case of the bound states. There exist two different possibilities 
depending on the sign of £ 2 . 



e 2 -M 2 

= rir 3 r 4 + r 2 r 3 r 4 + r 3 rir 2 + r 4 r\r 2 , 
S 2 

rir 2 r 3 r 4 . 



I e 2 — M l < , £ 2 = a(a 2 /M 2 ) < ; 



(2.5a) 



so, from (2.4) we conclude that two roots can be positive and two negative (or complex and 
conjugate to each other) 



n < , r 2 < , r 3 > , r 4 > . 



(2.56) 



that can be illustrated by Fig. 1. 

With the use of (2.4a) and (2.4b) one expresses the roots r\, r 2 through two classical turning 
points r 3 , r 4 

2ea T? 1 

H + r 2 = -77^ ^ - r 3 - r 4 , r 4 r 2 



n = 2 



7*2 = — 

2 



r 3 + r 4 



r 3 + r 4 



M 2 — e 2 
2ea 



e 2 — M 2 r 3 r 4 



M 2 -e 2 

2ea 
M 2 -e 2 



2ea \ 2 

7~3 + r 4 - TTo o I + 



4£ 2 



M 2 -e 2 J (M 2 -e 2 )r 3 r 4 



+ W h" 3 + r 4 - 



2ea 



4£ 2 



M 2 -e 2 / (M 2 -e 2 )r 3 r 4 



(2.5c) 



4 



To obtain two positive and two negative roots we require two conditions 

2ea 



— (ri + r 2 ) = r 3 + r 4 



2ea \ 2 
r3 + r A - — ) + 



M 2 -e 2 
4£ 2 



> 



> . 



M 2 - e 2 J ' (M 2 - e 2 )r 3 r A " " ' (2 ' M) 

The case with two positive and two complex roots (with two negative real parts) is realized if 
two following conditions are imposed 



-(n + r 2 ) = r 3 + r 4 - 



2ea 



M 2 



> 



r 3 + r 4 



2ea 



M 2 



+ 



4S 2 



(M 2 - e 2 )r 3 r 4 



< . 



(2.5e) 



->2. 



It should be noted that another possibility exists for the bound states at positive values for 



II 



e 2 — M 2 < , 



a 



from (2.4) we conclude that three roots can be positive and one negative 

ri < , r 2 > , r 3 > , r 4 > , 
r 2 r 3 r 4 , 2ea 



ri = 

r 3 r 4 + r 3 r 2 + r 4 r 2 

This case is illustrated in Fig. 2. 



, H = (r 2 + r 3 + r 4 
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< 
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(2.6a) 



(2.66) 



Figure 2: Finite classical motion: r G [0, r 2 ], r G [r 3 ,r 4 ] 
The formulas for two roots r\ and r 2 in terms of r 3 , r 4 are 
1 



ri 



?*2 = — 

2 



- r 3 + r 4 - 



2ea 



M 2 



2ea \ 2 
-\l[r 3 + ri- M2 _ e o 1 [ - 



4S 2 



(M 2 - e 2 )r 3 r 4 



r 3 + r 4 



2ea 



M 2 -e 2 



+ \/ [r 3 + r 4 



2ea 



+ 



4S 2 



M 2 -e 2 J (M 2 -e 2 )r 3 r 4 



but now the parameter is negative £ 2 > 0, and correspondingly the root r 2 is positive. 
The situation when one root is positive ant three are negative is also possible 



(2.6c) 



T\ < , r 2 > , r 3 < , r 4 > , 
ri r 2 r 3 



r 4 = -- 



that is illustrated in Fig. 3. 



r\r 2 + nr 3 + r 2 r 3 



P 2 (r) 



(2.6d) 



. r 4 



n r 2 r 3 



Figure 3: Finite classical motion: r G [0, r 4 ] 
Now let us consider possible infinite motions. The first possibility is 



III 



e 2 — M 2 > 



S 2 = o(a 2 /M 2 ) < ; 



n < , r 2 > , r 3 > , r 4 > . 
that is illustrated in Fig. 4. 

P 2 {r) 



n 



r 2 



r 3 •• r 4 



(2.7a) 
(2.76) 



Figure 4: Infinite classical motion: r G [r2,r 3 ],r G [r 4 ,+oo) 
Besides, another case can be realized if three roots are negative and one is positive 
ri < , r 2 < , r 3 < , r 4 > ; 
that is illustrated in Fig. 5. 



(2.7c) 



n r 2 r 3 

-• • • ■ 



P 2 (r) 



. r 4 



Figure 5: Infinite classical motion: r G [r4, +00) 

Finally, there exists one more case 

/y e 2 — M 2 > , £ 2 = fj(a 2 /M 2 ) > ; 
when two root are positive and two roots are negative 

n < , r 2 < , r 3 > , r 4 > ; 
that is illustrated in Fig. 6. 

P 2 (r) 



ri r 2 



7*3 ' 



r 4 



(2.8a) 
(2.86) 



Figure 6: Fig. 6. Infinite classical motion: R G [0,r2], r G [r4,+oo) 



3 Analytical treatment of the problem 

Let us turn back to eq. (2.2) and introduce a new variable 



i(-i^ 2 S 2 ) 1 / 4 r + S 
x = — — - ; — , r 



iT, (x + 1) 



(3.1a) 



i(-K 2 X 2 ) 1/4 r -S ' " (-K 2 S 2 )V4 {x- 1) • 

it is readily verified the main physical singularities r = and r = 00 in the new variable look as 

r = oo x = +l, r = x = -l. (3.16) 

In general, the variable x is complex-valued 

ivlr + £ 



x(r) = 



iAr - E ' 



( _ K 2 s2) l/4 = [(e 2 _ M 2 )S 2]l/4 = ^ 



(3.1c) 



There exist four different combinations of signs 

(e 2 - M 2 , S 2 ) =► J 



(+,+) 
(-,+) 



(+,-) 
(-,-) 



(3.2a) 



Below we will consider only the case of bound states when e 2 — M 2 < 0, which corresponds 
to variants (— , +) and (—,—). If X 2 (the case (—,+)), the variable x is complex- valued 



x(r) 



iAr + £ 



S 2 - ,4 2 r 2 



2AT,r 



iAr - S 
that is illustrated in Fig. 7. 



S 2 + A 2 r 2 £2 + ^2 r 2 



Im x 



Jf(r) 



x(r) 



(3.26) 



r = 00 Re x 



Figure 7: Variable x on complex plane 
If S 2 = — 6 2 < (the case (—,—)), the variable x is real- valued 

Ay _|_ h 

X ^ = A^b ' X e 1 _1 ' +1 ] • (3 - 2c) 
With the variable x, the differential equation (2.2) takes the form 

cte 2 + x-1 dx + \ ^/-K 2 Y?{x + l) 4 + y/-K 2 Y?(x - l) 4 ~ 
4ieaS ^_ 2ieaS - J^-^S 2 ) 1 / 4 

~ (-^2 S 2)l/4( x _ !)3 ~ ( x + 1)2 + (_K2 S 2)1/4( X _ X )2 + 

ieaS - J 2 (-^ 2 S 2 ) X / 4 -zeaS + J 2 (-if 2 E 2 )V 4 \ 

+ (_^ S 2)i/ 4 (x + i) + (-K 2 ^ 2 y/\ x - 1) J^- - (3 - 3) 

With the use of the notation 

(-i^S 2 ) 1 / 4 = A , -K 2 Y? = A 4 , A 2 = ^-K 2 Y? = ±iKZ , 
eq. (3.3) reads 

d 2 2 d ( 4A 2 _ AA 2 _4ieaE__ 

d!?* + + ^"(x + l) 4 ~ (x - l) 4 " A(x - l) 3 ~ 
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J 2 2ieoiS - J 2 A iea^-J 2 A -ieaL + J 2 A\ 
~(x + l) 2+ A(x-l) 2 + A(x + l) + A(x-l) )^ x >-°- & A > 

Applying the substitution 

f(x) = (x + l) B (x - if exp + ^) (3.5) 

we arrive at the equation for F 

d 2 F ( 2B 2C + 2 D D \ df 



dx 2+ \x + l + x-l (X + 1) 2 (x-l) 2 ) dx + 
'D 2 -16A 2 D 2 -16A 2 



+ 



4(x + l) 4 4(x-l) 4 



D{B - 1) CT> ,4 + 4ieaE D 2 + 8£ 2 - 85 + 4CL> + 4D - 8 J 2 
~ (x + 1) 3 A(s - l) 3 + 8 (x + 1) 2 + 

D 2 A + 8(L4 - 4.BD A + 8C 2 A + lUeaE - 8 J 2 A 
+ 8A(x-l) 2 + 

D 2 A + 2CD A - 2BD A - 8BC A + 2DA-8BA + 8ieaT, - 8 J 2 A 
+ 8A{x + 1) + 

-D 2 A - 2CD A + 2BD A + 8BC A-2DA + 8BA- 8ieaS + 8 J 2 A\ 



8A(x - 1) 



F = 0. 



When B, C, D are given as (see (3.5)) 



1 „ 1 AiKY, 

B = - C = — , D = ± = ±4A , 

2 2 ' ,4 



eq. (3.6) becomes more simple 
d 2 F 



dx 2 



x-l *\(x + l)(x-l) 



1 1 D D 

+ 



+: 



.x + 1 x-l (x + 1) 2 (x-l) 2 

1 



dF 
dx 



(D 2 A - 8J 2 A -2A- 16iea£)x 2 + 



(3.6) 



/ = J P\ exp f ) F(x) (3.7) 



2,4(2; + 1) 3 (x-l) 3 

+ {8DA - 32iea£)x + 2A- 16ieaS - L> 2 ,4 + 8J 2 vl 1 F = . (3.8) 

It coincides with the double confluent Heun equation [3j|4] for -ff(/i, /3, 7, 5, z): 

d 2 H f 1 1 n /I \dH 

dx 1 " + \x + l + x - 1 ~ 2(x + l) 2 ~ 2(x - 1) 2 J ~dx~ + 
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with parameters 



^±^±M^±l H = (3 . 9o) 

WieaY, 



^ = 2D = ±8 A , 7 



A 



f3 = -1 - 4J 2 + 8A 2 - 8iea^- , 

A 

5 = +1 + 4J 2 - 8A 2 - 8iea- ; (3.96) 



with additional constrain (3 + S = 7 . 



4 Numerical simulations 

To possibility for quantum mechanical bound states with energy ei[— M, +M] there must corre- 
sponds in classical description a finite region for classical motion. Let us examine the condition 
at which the 4-nd order polynomial (this analysis will enable us to find left limiting boundary 
V for possible quantum energy levels, e > eo): 

m , _ (e 2 ~ M 2 )r 4 + 2ear 3 - j(J + l)r 2 + aa 2 /M 2 

U V r J £ 2 _ ^2 U { ' 

just starts to have a double root ro at positive real axis 

II(r) = (r - r ) 2 (r - a + i&)(r - a - 16) = . (4.2) 

This is just the bifurcation value of e. 

It is convenient to measure e in unit of M introducing e = eM. From comparing (4.1) and 
(4.2) we get the system of algebraic equations for ro, a, b, e 

- (a 2 + iV 2 + _i^__ = , 

2r (a 2 + b 2 + ar ) = , 



From whence it follows 



M(ro — e 2 ro) — ea 



(-l + e 2 )M 

2 _ (ea((-l + e 2 )Mr + 

' ~ (-l + e 2 ) 2 M 2 ' 

3ea + v/8(-l + e 2 )j(l + j) + 9e 2 a 2 
r ° ~ 4(-l + e 2 )M 
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and we arrive at a rather complicated equation for e 

— 16(— 1 + e 2 ff - 



l + e 2)2 j 4_ 27e 4 a 4. 



9e 3 a 3 



-1 + e 2 )j(l + j) + 9e 2 a 2 - 



2ea 



4e(-l + e 2 )ja(9ea + 2 v /8(-l + e 2 )j(l + i) + 9e 2 a 2 )- 
4(-l + e 2 )j 2 (-2 + e 2 (2 + 9a 2 )+ 
-1 + e 2 )j(l + j) + 9e 2 a 2 ) + 32(-l + e 2 ) 3 a 2 a = 



There exists additional constraint (condition for real-valuedness of ro( 



e > 



8j + 



8j + 8j 2 + 9a 2 



(4.4) 



(4.5) 



Numerically one easily finds the value of e m j n which is a lower boundary for the existence 
of bound states 

At j = 0, cr = — l,a = 1,M = 1 we obtain e m j n = 0.614659, the plot of curve is illustrated 
in Fig. 8. 



15 



10 - 



Figure 8: U(r) at j = 0, <r = -1, a = 1, M = 1 

As one can find by numerical simulations at a > possible values for e lay outside the interval 
[-1,1] and lead to ro < 0. 

After obtaining information of regions of possible bound state energy e for the system with 
small quantum number j, . . . we can start to construct numerical solution for the lowest eigen- 
state. To this goal we extend the known Ritz variational approach [5J, performing simulations 
for simplicity at fixed values M = l,a = — 1 , j = . 



11 



First we find the asymptotic behavior of the eq. (1.7) at the origin. Selecting the most 
singular terms we get 

f"(r) - a 2 f(r)/r 4 = 

with the leading term in asymptotic solution of the form r<l, f(r) = l/(2a) exp(— a/r)r. In 
a similar way at infinity we introduce the variable u = 1/r, rewriting the equation and selecting 
most singular terms only we get 

y"{u)-y{u)/u 4 (l-e 2 ) = (4.6a) 

with leading term in asymptotic solution at r — > oo as 

f{r) ~ r- 1 (exp(-(v / l - e 2 r)/(2y/l - e 2 ). (4.66) 

Then the trial lowest eigenstate can be chosen as a product of these leading terms with yet 
unknown energy. Introducing notion k = y/l — e 2 first we find the normalization condition 

J (exp(- (a/r) - r «) / (4«)) 2 r 2 dr = ^ ^ffi^ (4.7) 

where Ki is the Bessel function of imaginary argument of order i. 

Multiplying from the left the left-hand side of eq.Q on normalized trial function f(r) and 
integrate over radial coordinate. As a result we get the expression of the form 

K 3 (4^ok) (an - k 2 + e 2 - l) + a(2e - 1)K 2 (4 v / o^) = (4.8) 

which is a second order algebraic equation on e with elementary solution for roots. 

Now, in accord with usual Ritz variational method we have to minimize one of the root 
on respect to k that leads to the approximate value of the lowerst eigenstate energy. The 
appropriate dependence is shown in Fig. 9 for a = 1. 




Figure 9: The dependencies of the roots £i(k), e 2 (k) 



As one can see from Fig. 9 we indeed have a minimum for the second root, its value and 
location are e = 0.749279 at k = 0.625342. The appropriate eigenfunction and squared redial 
momentum P 2 (r) are shown in Fig. 10 
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Figure 10: Plots of the lowest energy eigenfunction and P 2 (r) at M = 1, a = = 0, a = 1 
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